Glauber-model analysis of total reaction cross sections for Ne, Mg, Si, and S isotopes 

with Skyrme-Hartree-Fock densities 
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A systematic analysis is made on the total reaction cross sections for Ne, Mg, Si, and S isotopes. 
The high-energy nucleus-nucleus collision is described based on the Glauber model. Using the 
Skyrme-Hartree-Fock method in the three-dimensional grid-space representation, we determine the 
nuclear density distribution for a wide range of nuclei self-consistently without assuming any spatial 
symmetry. The calculated total reaction cross sections consistently agree with the recent cross 
section data on Ne-|-^^C collision at 240^4 MeV, which makes it possible to discuss the radius and 
deformation of the isotopes. The total reaction cross sections for Mg-f ^^C, Si-|-^^C and S-|-^^C cases 
are predicted for future measurements. We also find that the high-energy cross section data for O, 
Ne, and Mg isotopes on a ^'^C target at around 1000 AMeV can not be reproduced consistently with 
the corresponding data at 240 AMeV. 
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I. INTRODUCTION 

Advances in measurements of unstable nuclei have 
been providing information on exotic nuclei toward the 
neutron and proton driplines. Recently total reaction or 
interaction cross sections have been measured in p/-shell 
region [ll-Q . They exhibit interesting trends which imply 
exotic structure, for example, halos, skins, deformations 
etc., as we approach the driplines. In particular, the first 
observation of a halo structure of '^^Ne [1] has stimulated 
several theoretical studies 0-0| • 

Nuclear deformation is one of the unique properties of 
a finite quantum system. If the intrinsic wave function 
shows some deformation, the nuclear radius becomes ef- 
fectively large because the ground-state wave function is 
expressed as a superposition of the intrinsic wave func- 
tions with different orientations. Since the total and in- 
teraction cross sections are closely related to the nuclear 
radius, it is interesting to investigate a relation between 
these cross sections and nuclear deformation. 

In this paper, we present a systematic analysis of to- 
tal reaction cross sections, on a ^^C target, of unstable 
nuclei in the sd- and p/-shell region, focusing on nuclear 
size properties, especially radius and deformation, and 
their relationship with the cross section. This study is 
motivated by the recent systematic measurements of the 
total reaction cross sections of the p/-shell region [ll-Q . 
We also predict the cross sections of other nuclei for fu- 
ture measurements. The analysis will give us important 
information on the interplay between the nuclear struc- 
ture and the cross section. 

In Refs. [1, @, two of the present authors (W.H. and 
Y.S.), B. Abu-Ibrahim, and A. Kohama performed sys- 
tematic analyses of the total reaction cross sections of C 
and O isotopes on a ^^C target. In these studies, the 
wave functions generated from a spherical Woods-Saxon 
model were employed and the center-of-mass motion was 



appropriately removed. The Glauber model, which gives 
a fair description of high-energy nucleus-nucleus colli- 
sions from nucleon degrees of freedom, provided a good 
agreement with experiment. However, in the present pa- 
per, we study the region including the island of inver- 
sion where the nuclear deformation may play a significant 
role. Thus, we perform the Skyrme-Hartree-Fock calcula- 
tion in the three-dimensional (3D) coordinate-mesh rep- 
resentation which allows us to treat the correct asymp- 
totic behaviors of single-particle orbitals and the nuclear 
deformation self-consistently. Although the parameter 
sets of Skyrme energy functionals are mostly determined 
by fitting the properties of heavy closed-shell nuclei and 
the nuclear matter, they are known to give a good de- 
scription of light nuclei as well. This universality is one 
of the main advantages of the density-dependent mean- 
field models. In order to see the sensitivity of the nuclear 
deformation to the radius, two different types of Skyrme 
parameter sets are employed for comparison. 



This paper is organized as follows. In Sec. [Hi we briefly 
explain our reaction and structure models to calculate 
the total reaction cross sections. Sec. HH] presents the 
calculated total reaction cross sections on a -^^C target 
at the medium energy of 240 AMeV, corresponding to 
the recent experiments carried out at RIKEN. We com- 
pare our theoretical cross sections and the experimental 
ones for Ne in Sec. IIII Al Motivated by a very recent 
measurement we show our results for Mg isotopes in 
Sec. IIII B1 Predictions for Si and S isotopes are also made 
in Sec. IIII Cl We discuss in Sec. IIII Dl a problem concern- 
ing the total reaction cross section data at around 1000 
AMeV for O, Ne, and Mg isotopes. Conclusion is drawn 
in Sec. HVl 
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II. MODELS 
A. Glauber model 

We describe a high-energy nucleus-nucleus collision in 
the Glauber formalism • The total reaction cross sec- 
tion is obtained by integrating the reaction probability 
P{b) over the impact parameter b 



OR = j dbP{b), 



with 



P{b) = 1 



(1) 



(2) 



where is the phase-shift function for the elastic scat- 
tering of a projectile (P) and a target (T). The phase 
shift function is given in terms of the ground-state wave 
functions of the projectile and the target as 
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where Si is the transverse component of the ith nucleon 
coordinate and Tnn is the nucleon-nucleon profile func- 
tion which we parametrize in the following form 
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and 13 NN used here are 



The parameter set of cr* 
given in Ref . [ll| . 

Though the multiple integration in Eq. ([Sj may be per- 
formed with a Monte Carlo integration as was done for 
light nuclei [T^], wc take a usual approximation in this 
paper. The optical limit approximation (OLA) offers the 
most simple expression, which is obtained by taking only 
the first-order term of the cumulant expansion [l^ of Eq. 
©, that is 
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which involves a double-folding procedure of the projec- 
tile and target densities {pp, pr) with the effective inter- 
action r^rTv. The OLA misses some higher-order terms 
of F^vAf and multiple scattering effects. The resulting to- 
tal reaction cross section of the nucleus-nucleus collision, 
tends to overestimate the measured cross section 

In Ref. [ll] , another expression of evaluating the phase 
shift function is proposed in order to take account of 
the multiple scattering processes missing in the OLA. To 
derive the formula, first we use the cumulant expansion 
of the phase shift function for the nuclcon-targct system 



and take the first-order term back to the original phase 
shift function ([3]). Using the cumulant expansion again, 
its leading term is called the Nucleon- Target formalism 
in the Glauber model (NTG) approximation, which is 
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In numerical calculations shown in Sec. IIIIl we use the 
symmetrized expression of Eq. ([6]) by exchanging the 
role of the projectile and target nuclei. Note that the 
NTG approximation requires the same inputs as those of 
the OLA. It gives a simple but fair description for high- 
energy reactions. For example, the total reaction cross 
sections of ^^C+^^C collisions are improved very much in 
a wide energy range Q. 



B. Skyrme-Hartree-Fock method in 3D 
coordinate-mesh representation 

We perform the Skyrme-Hartree-Fock calculation for 
the density distribution of a variety of projectiles. The 
ground state is obtained by minimizing the following en- 
ergy density functional 



E[p] = E. 



N 



En — Er 
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For the ground states of even-even nuclei, the nuclear 
energy E^ is given by a functional of the nucleon den- 
sity pq{r), the kinetic density Tq{r), the spin-orbit- 
current density V ■ Jq{r) {q = n,p). The Coulomb 
energy Ec among protons is a sum of direct and ex- 
change parts. The exchange part is approximated by 
means of the Slater approximation, cx J drpp{r)^^^ . The 
center-of-mass recoil effect is also treated approximately 
by a subtraction of the expectation value of £^cm = 
(E.pOV(2mA)«E,:P?/(2™A). ^ 

Every single-particle wave function (j)i{r,a,q) is rep- 
resented in the 3D grid points with the mesh size of 
Ax = Aj/ = Az = 0.8 fm. All the grid points inside the 
sphere of radius of 15 fm are adopted in the model space. 
The ground state is constructed by the imaginary-time 
method [TBj with the constraints on the centcr-of-mass 
and the principal axis: 



j drxp{r) = j dryp{r) = j drzp{r) = 0, 
j drxyp{r) — J dryzp{r) — J drzxp{r) = 



(8) 



For odd-A nuclei, we adopt the half-filling approxima- 
tion, thus the time-odd densities, such as the current den- 
sity and the spin density, do not contribute to the present 
calculation. The computer program employed in the 
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present work has been developed previously for linear- 
response calculations [iGl - flsj . including all the time-odd 
densities. However, those time-odd parts automatically 
vanish in converged self-consistent solutions. Because of 
this simple filling treatment for odd-A nuclei, we do not 
expect a quantitative description of odd-even effect in 
our calculation. For nuclei very near the neutron drip 
line, the nuclear radius becomes extremely sensitive to 
the neutron separation energy. In other words, a light 
modification of the Skyrme parameter set could drasti- 
cally modify the final result. Thus, we do not show such 
results for the drip-line nuclei in Sec. IIIIl 

We calculate isotopes with even proton numbers, Z = 
8jl0, 12, 14, and 16. The two parameter sets of SkM* 
[l9l | and SLy4 [l^] are employed in the following, and we 
compare results obtained with these two parameter sets. 
The SkM* functional is known to well account for the 
properties of the nuclear deformation, while the SLy4 is 
superior to SkM* in reproducing the total binding energy. 

The self-consistent solution becomes, in most cases, a 
deformed intrinsic state \^k) with the definite K quan- 
tum number. The ground-state wave function in the 
laboratory frame is constructed according to the strong- 
coupling scheme [2l| as 
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+ (-l)'+^>A'X-A'H}, 

where DI^^{uj) is the D function depending on the three 
Euler angles w. The density distribution in the ground 
state with / = _ftr is simply given by 



(10) 



taking the average over states with different magnetic 
quantum number M. The density p"^'^'°\r) of Eq. ([TU]) 
is utilized in Eq. (|6]) as the projectile density distribution 
ppir). Since the Skyrme energy functional is constructed 
so as to give a density distribution p{r) in the center-of- 
mass frame, we directly use the density in Eq. (|10p in 
the NTG formula ^. 



A. Neon isotopes 

Figure [1] displays the point matter, neutron, and pro- 
ton root- mean-square (rms) radii of Ne isotopes as a func- 
tion of neutron number N . The enhancement of the mat- 
ter radius in the neutron-rich isotopes is dominantly due 
to that of the neutron radius. An interesting observation 
is that both the neutron and proton radii follow the same 
behavior as the matter radius. Both SkM* and SLy4 en- 
ergy functionals produce similar results for N < 18. We 
see, however, noticeable difference beyond TV = 18, and 
then two results coincide again at TV = 24. What kind of 
structure property is responsible for these differences? 
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FIG. 1: (Color online) Point matter, neutron, and proton 
rms radii of Ne isotopes calculated with the SkM* and SLy4 
interactions. 



III. RESULTS AND DISCUSSION 

Here we treat a reaction on a ^^C target. Proton and 
neutron density distributions of ^^C are obtained in the 
same way as in Ref. . The obtained density reproduces 
the proton charge radius determined by an electron scat- 
tering. The total reaction cross sections of both p+^^C 
and ^^C-|-^^C are reproduced very well in a wide energy 
range d, [ll| . We adopt this density distribution for the 
target nucleus. The density distribution of the projec- 
tile nucleus is determined by the Skyrme-Hartree-Fock 
calculation as explained in Sec. IIIBI For both the pro- 
jectile and target, we use proton and neutron densities 
separately. 



Figure [21 shows the calculated neutron Fermi energies. 
Both SkM* and SLy4 results exhibit similar behavior, 
except for very neutron-rich isotopes with N > 20, in 
which the neutron separation energy is calculated to be 
larger for SkM* than SLy4. Since a smaller separation 
energy is expected to produces a larger matter radius, the 
behavior of the Fermi energy or the neutron separation 
energy does not explain the difference of the calculated 
matter radii for nuclei with 18 < < 24 (Fig. [J). On 
the contrary, this gives an opposite effect. The matter 
radius is predicted to be larger for the SkM* calculation. 

The different matter radii arc found to be mainly due 
to the difference in the nuclear deformation. We calculate 
the quadrupolc deformation parameters from the mass 
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FIG. 2: (Color online) Neutron Fermi energies of Ne isotopes 
as a function of neutron number A^. 



quadrupole moments 



/5 = \//3|o + /3|2, 



(11) 



with 



/320 
/322 



n{2z^ - - 
5 (r2) 



(12) 



The quantization axis is chosen as the largest (smallest) 
principal axis for prolate (oblate) deformation. Figure [3] 
presents the quadrupole deformation parameters /3 ob- 
tained from the Hartree-Fock (HF) solutions. The pos- 
itive (negative) values of /3 indicate the prolate (oblate) 
deformations. Beyond N = 16, we find a significant dif- 
ference in the magnitude of /?, between SkM* and SLy4. 
In the SkM* calculation, the deformation parameter (3 
is positive and rapidly increases for N > 18. On the 
other hand, the SLy4 predicts the nearly spherical ground 
states (weakly oblate) for TV < 20 and changes its shape 
into prolate for N > 20. The value of \/3\ is always larger 
in SkM* for 18 < iV < 24. This nicely corresponds to 
the observation of the different behavior in the matter 
radii at 18 < < 24 in Fig. [TJ The matter radii co- 
incide again at = 24, so does the predicted /3 values. 
It should be noted that the kink behavior in the matter 
radii at = 20 (SLy4) in Fig. [T] can be explained by the 
onset of the deformation at iV > 20, while the kink at 
= 14 is due to the occupation of the Si/2 orbitals at 
N > 14. 

The proton number Z = 10 of Nc isotope is known to 
strongly favor the prolate shape, because it corresponds 
to the magic number at the superdeformation [la . [2l| . 



This is consistent with the results in Fig. [3] for Ne iso- 
topes with A^ ~ 10. Recent experiments reveal that the 
excitation energy of the first J'^ ~ 2+ states E{2f)_m 
even Ne isotopes decrease from Af = 16toA^ = 22 [23 |. 
This is qualitatively consistent with the increase of the 
quadrupole deformation calculated with SkM*. 
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FIG. 3: (Color online) Quadrupole deformation parameter (5 
of Ne isotopes calculated with the SkM* and SLy4 interac- 
tions. 
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FIG. 4: (Color online) Total reaction cross sections of Ne 
isotopes on a ^^C target at 240 ^MeV. Experimental data 
are the interaction cross sections at 240 AMeV taken from 
Ref. 0]. 



The calculated total reaction cross sections of Ne iso- 
topes incident on a "'^^C target at 240 A MeV are shown 
in Fig. 21 as a function of neutron number A^. The agree- 
ment with the recent data is very good. The slope of 
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the curve in Fig. |4] changes at A'' = 14. The cross section 
increases gradually from N = 10 to 14 and shows a rapid 
rise from = 14. This neutron-number dependence of 
the total reaction cross sections also well corresponds to 
the calculated matter radii. 

For the nuclei with = 10 14, the calculated 
cross sections are slightly larger than the observed values, 
while the agreement between the theory and experiment 
is fairly good in the neutron-rich region beyond N ~ 14. 
This may be explained by the difference between the total 
reaction and interaction cross sections. The theory gives 
the total reaction cross section cr/j while the experiment is 
the interaction cross section cr/ . Since the latter cross sec- 
tion does not contain inelastic reaction cross sections, the 
inequality > aj always holds. Contributions of the 
inelastic processes are classified as two cases: (i) The pro- 
jectile nucleus is excited to particle-bound excited states 
that are located below particle-emission threshold, while 
the target can take any states, (ii) The target nucleus is 
excited, but the projectile remains in the ground state. 
The difference between the total reaction and interaction 
cross sections is expected to be smaller as going closer to 
the neutron drip line. The difference is estimated to be 
at most 100 mb in a phenomenological way (23| . 

For the neutron-rich nuclei ^^~^^Ne, the results of 
SkM* nicely reproduce an average behavior of the exper- 
imental cross sections, but not their odd-even staggering. 
This may be due to the pairing correlation [23|, which is 
neglected in the present study. Our result for these nuclei 
is similar to the one obtained with the deformed Woods- 
Saxon potential [25|. Approaching the neutron dripline, 
the neutron separation energy becomes very small. In 
such a situation, a spatial extension of the last neutron 
orbit is very sensitive to the separation energy. The 3D 
coordinate-space representation for the single-particle or- 
bitals is suitable for the treatment of their asymptotics. 
However, in the case of '^^Ne, the measured neutron sep- 
aration energy is only about 0.33 MeV with a large error 
bar [2^ , while the calculated separation energy is as large 
as 2 MeV with the SkM* interactions. Together with the 
half-filling approximation, the present Skyrme function- 
al may not provide us with such a precise description of 
a loosely bound system. To artificially adjust the Fermi 
level to 300 keV, we multiply the mean-field potentials 
{Un and Up [l^) obtained with SkM* by a common scal- 
ing factor, then, the total reaction cross section is easily 
increased by 140 mb. This, in fact, significantly overesti- 
mates the experimental cross section. Thus, an accurate 
measurement of the neutron separation energy is highly 
desired. 

Figurc[5]presents the calculated density profiles of ^"Nc 
in both linear and logarithmic scales. The nucleus ^°Ne 
is strongly deformed with SkM* (/? ~ 0.35), while the 
shape is almost spherical with SLy4 (/? = —0.06). In this 
case, two parameter sets exhibit quite different density 
profiles, especially in the outer region that is crucial to 
determine the matter radius. The density calculated with 
SkM* is larger than the one with SLy4 at r > 5 fm. 



0.2 



0.15 



0.05 





1 1 


1 1 
SkM* 






CI wA 

Ol— yr 




v 

v 

V 
v 
V 






1 1 






FIG. 5: (Color online) Upper panel: Point matter density 
distributions of "^"Ne calculated with the SkM* and SLy4 in- 
teractions. Lower panel: The same as the upper panel but 
drawn in a logarithmic scale. 



In order to see how the above density difference leads 
to the difference in the total reaction cross section, we 
plot the difference of the reaction probabilities calculated 
with the SkM* and SLy4 interactions 



(13) 



Note that D{b) depends only on 6 = In Fig. [SI we 
show D{b) multiplied by 27r6 for '^°Ne projectile on both 
^^C and proton targets for comparison. Though we see 
some difference in the interior density as shown in Fig. [SJ 
it has almost no contribution to D{b) for a ^^C target, 
while D{b) on a proton target shows oscillatory behavior 
reflecting the difference of the density profiles between 
2^4 fm. The reaction on a ^^C target is insensitive to 
the internal density profiles. In contrast, we can see the 
difference in the outer region (Fig. The reaction on a 
^^C target is advantageous to probe the density at large 
distances. 



B. Magnesium isotopes 

Let us next discuss Mg isotopes. Figures [7] and [5] 
display the quadrupole deformation parameters and the 
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FIG. 6: (Color online) Impact parameter dependence of the 
difference of the reaction probabilities for '^"Ne calculated with 
the SkM* and SLy4 interactions. 
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matter, neutron, and proton rms radii for Mg isotopes. 
The trend of the deformation parameters is very similar 
to that of Ne isotopes. Similarly to Ne case, the /? val- 
ues obtained with SkM* are larger than that with SLy4 
for 18 < < 24. This behavior influences the matter 
radii as well as the neutron and proton radii as shown in 
Fig. [S] We also plot the radii obtained by the fermionic 
molecular dynamics approach (FMD) [2^ [2^. Though 
the FMD radii tend to be slightly smaller than our HF 
radii in 12 < iV < 18, the HF radii with SkM* are very 
close to the FMD radii for N > 18. As shown in Ref. Hj, 
the trend of the quadrupole deformations is also similar 
to that obtained with SkM*. 

The deformation /3 calculated with SkM* is very sim- 
ilar to that of Ne case. The difference can be observed 
only for iV = 14 ~ 16: the Ne isotopes have oblate shape 
with small deformation, while the Mg isotopes have pro- 
late shape. Experimental excitation energies of the first 
2+ states, E{2+), of ^s.soj^g ^ ^^g^ ^8) arc almost 
the same, whereas we see a rapid decrease of i?(2^) in 
26,28]\jg |22[. Again, this is qualitatively consistent with 
the calculation. 

In Fig. [3] we show the total reaction cross sections cal- 
culated using the SkM* and SLy4 interactions for Mg 
isotopes incident on a ^^C target. A comparison with ex- 
periment is very interesting to judge whether the neutron 
rich Mg isotopes favor the strong deformation or not and 
also whether the HF densities at 12 < < 18 can re- 
produce experimental cross sections or not because they 
predict larger radii than the FMD calculation. Here we 
comment on the cross section of ^^Mg. Since the neu- 
tron separation energy of '^^Mg is considered to be small 
(-162 keV) [131, and the /? value of ^^Mg is fairly large, 
a situation similar to ^^Ne may occur in "^^Mg as well. 
If that is the case, the cross section predicted for "^^Mg 
will be increased than the present value. We hope the 
measured values of the total reaction (interaction) cross 



FIG. 7: (Color online) Quadrupole deformation parameter /3 
of Mg isotopes calculated with the SkM* and SLy4 interac- 
tions. 



C. Silicon and sulfur isotopes 

We now discuss the cases of Si and S isotopes, which 
will show deformation different from Ne and Mg isotopes. 
Figures [TU] and [TT] plot the quadrupole deformation pa- 
rameters as well as the point matter, neutron, and proton 
rms radii for Si and S isotopes, respectively. We also see 
some kinks in the calculated matter, neutron, and pro- 
ton radii in Fig. [TTJ corresponding to the peaks of the 
quadruple deformation parameter displayed in Fig. [TOl 
The magnitudes of /? for Si and S isotopes are signifi- 
cantly smaller than those of Ne and Mg isotopes, in the 
region of iV w 20, which suggests that the island of in- 
version does not reach Z = lA (Si). This is probably 
due to that 14 protons favor the oblate shape, while the 
neutrons tend to deform the system into a prolate shape 
for N > 16. Thus, the competition between protons and 
neutrons determines the deformation. The enhancement 
of the matter radii for > 14 is not as strong as that 
in Ne and Mg isotopes, corresponding to the moderate 
change of /3 values. The matter radii of Si and S iso- 
topes show similar enhancement. The proton radii of S 
isotopes, however, tend to be larger than those of Si iso- 
topes due to the occupation of the s orbit. The neutron 
and matter radii of S isotopes also show kink behaviors 
at = 28. This may be due to the decrease of the de- 
formation at A^ = 26 to 28 and the occupation of \ow-i 
{p) orbitals beyond TV = 28. Note that similar kink be- 
haviors at A" = 28 are predicted in the low-energy El 
strength for S isotopes but not for Si isotopes [l^ . 

In Fig. [1^ we plot the predicted total reaction cross 
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FIG. 8: (Color online) Point matter, neutron, and proton 
rms radii of Mg isotopes calculated with the SkM* and SLy4 
interactions. The radii obtained by the FMD [2^, [2^ are also 
shown. 
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sections of Ne and Mg isotopes. In fact, the increase of 
ur from iV = 18 to 24 is 210, 180, 140, and 130 mb, for 
Ne, Mg, Si, and S isotopes, respectively. 
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FIG. 10: (Color online) Quadrupole deformation parameters 
of Si and S isotopes calculated with the SkM* interaction as 
a function of neutron number A^. 
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FIG. 11: (Color online) Point matter, neutron, and proton 
rms radii of Si (left) and S (right) isotopes as a function of 
neutron number A'^. The SkM* interaction is used. 



FIG. 9: (Color online) Total reaction cross sections of Mg 
isotopes incident on a ^^C target at 240 AMeV. 



D. Problems in the high-energy data 



sections of Si and S isotopes on a ^"^G target at 240 A 
MeV. For the sake of comparison, the ur values for Ne 
and Mg isotopes are also drawn. Though some behavior 
indicating the nuclear deformation appears to persist in 
Si and S isotopes as well, the change of the cross sections 
in 18 < iV < 24 is not as drastic as that of the cross 



We have shown that the present theory reproduces 
the total reaction or interaction cross section data at 
the medium energy of 240v4MeV. Since it is based on 
the Glauber model, our theory should provide us with 
a very good description of high-energy collisions around 
1000 AMeV. In fact, as shown in Ref. Q, the total reac- 
tion cross section of ^^C-I-^^C collision is reproduced very 
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FIG. 13: (Color online) Total reaction cross sections of O 
isotopes on a ^^C target at 1000 AMeV. The SkM* interac- 
tion is used. Measured interaction cross sections at around 
950 AMeV are taken from Ref. [s^l (open diamonds). Most 
recently measured interaction cross sections by Ref. [s^ are 
denoted by filled diamonds. 



FIG. 12: (Color online) Total reaction cross sections of Ne, 
Mg, Si and S isotopes on a ^^C target at 240 AMeV as a 
function of neutron number A*'. The SkM* interaction is used. 
The cross sections for Ne, Mg, and Si isotopes are plotted by 
subtracting 300, 200, and 100 mb, respectively. 



well in the wide energy ranging up to 1000 AMeV: The 
calculated cross sections at 240 and 1000 AMcV are 790 
and 850 mb, respectively, which is in excellent agreement 
with the measured total reaction (interaction) cross sec- 
tions, 782±10mb [HI at 250 AMeV and 853±6mb [H 
at 950 AMeV. We point out in this subsection, however, 
that some cross sections for O, Ne, and Mg isotopes at 
around 1000 AMcV are difScult to reproduce. 

We start with a problem concerning the interaction 
cross sections of O isotopes. The anomalously large cross 
section of ^^O on a ^^C target at around lOOOAMeV has 
been a long standing problem js^ - ls^ ]. A large jump of 
the measured cross section tj/ from ^^O to ■^'^O was an- 
alyzed by a model of adding one neutron in a loosely 
bound Is orbit to the ^^O core [s^. However, it failed 
to explain the data because one neutron separation en- 
ergy (2.74 MeV) is too large to form a halo structure. 
Figure [T3] displays the total reaction cross sections of 
O isotopes calculated with the SkM* interaction. Mea- 
sured interaction cross sections are taken from Ref. [lH 
for 13-240 and from Ref. [s^ for 22,230. The latter are 
new data remeasured in 2011. In the case of ^^O, our 
result agrees well with the new data [s^ as well as the 
calculation based on a phenomenological Woods-Saxon 
potential @. The former experimental value for 2^0 (32j 
seems to be too large. 

Next, we show in Fig. [T3]the total reaction cross sec- 



tions of Ne and Mg isotopes on a -^^C target at around 
1000 AMeV. Experimental cross sections are taken from 
Refs. [l], lllj. Again, a considerable difference between 
the theory and experiment is observed. This difference 
is unexpected, however, in view of the fact that we have 
a good agreement with the measurement at 240 AMeV 
for Ne isotopes (Fig. |4]). It is also noted that, although 
the calculation produces a smooth neutron number de- 
pendence, the measured cross sections show an irregular 
decrease at ^^Mg. 

Since our calculations at 240 AMeV reproduce the re- 
cent measurement very well and since the total reaction 
cross sections are expected to be very close to the interac- 
tion cross sections at around 1000 AMcV, we think that 
updating experimental cross section data is necessary to 
resolve the above problems. 



IV. CONCLUSION 

We have made a systematic analysis of the total reac- 
tion cross sections of Ne, Mg, Si, and S isotopes on a ^^C 
target. The densities are obtained using the Skyrme- 
Hartrcc-Fock method on a full three-dimensional grid 
space. High-energy nucleus-nucleus collisions arc de- 
scribed using the Glauber model. Both structure and 
reaction models employed here have no adjustable pa- 
rameter. 

Comparing our results with the recent measurements 
of the total reaction cross sections for Ne isotopes at 
240 AMeV, we find that a good agreement is obtained 
for both stable and unstable projectiles. We have shown 
that the nuclear deformation plays an important role to 
determine the matter radius of the neutron-rich isotopes. 
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The systematic measurements of an (tj/) for a long chain 
of isotopes may reveal the enhancement of the nuclear 
size which can be a signature of the nuclear deformation. 
We have found that the proton radii also follow the same 
behavior of the matter radii. If one can measure the 
proton radii, for example, by charge changing cross sec- 
tions, they also give information on nuclear deformations 
as well as nuclear skin-thicknesses. 

We have predicted the total reaction cross sections for 
Mg, Si and S isotopes on a ^^C target. We have also 
pointed out some contradictions in the high-energy cross 
section data at around 1000 ^MeV. In spite of the fact 
that we can excellently reproduce the experimental data 
at 240 ^MeV, we observe considerable disagreement be- 
tween our calculation and the measured cross sections 
for O, Ne, and Mg isotopes on a ^^C target at around 
1000 AMeV. Resolution of these problems requires fur- 
ther investigations in both theory and experiment. 

It is interesting to examine the total reaction cross sec- 
tions on a different target because they give different sen- 
sitivity to the density profile of the projectile nucleus. 
Since a ^^C target has a finite size, the reaction mainly 
takes place at the nuclear surface where more particles 
are involved. On the contrary, a proton target, as a point 
particle, is able to probe the inner region of the nucleus. 
This will provide us with information on the interior den- 
sity differently from the carbon target. The work along 
this direction is underway. 



FIG. 14: (Color online) Upper panel: Total reaction cross 
sections of Ne isotopes on a ^'^C target at lOOOylMeV. The 
SkM* interaction is used. Measured interaction cross sections 
at around 950AMeV are taken from Ref. [32|. Lower panel: 
Total reaction cross sections of Mg isotopes on a carbon tar- 
get at gOOAMeV. The SkM* interaction is used. Measured 
interaction cross sections at around 950 AMeV are taken from 
Refs. [s^l (closed diamonds) and [J] (filled diamonds). 



A similar trend in the total reaction cross sections is ex- 
pected for Mg isotopes whose deformation behaves like 
that of Ne isotopes. We have also shown that the to- 
tal reaction cross sections on a ^^C target is sensitive to 
the nuclear density, especially near its surface where the 
deformation contributes to changing the density profile. 
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